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ABSTRACT 

We investigate the instability of purely poloidal magnetic fields in nonrotating neutron stars by means of 
three-dimensional general-relativistic magnetohydrodynamics simulations, extending the work presented in 
ICiolfi et all d201 ll) . Our aim is to draw a clear picture of the dynamics associated with the instability and to 
study the final configuration reached by the system, thus obtaining indications on possible equilibria in a mag- 
netized neutron star. Furthermore, since the internal rearrangement of magnetic fields is a highly dynamical 
process, which has been suggested to be behind magnetar giant flares, our simulations can provide a realis- 
tic estimate of the electromagnetic and gravitational-wave emission which should accompany the flare event. 
Our main findings are the following: (i) the initial development of the instability meets all the expectations 
of perturbative studies in terms of the location of the seed of the instability, the timescale for its growth and 
the generation of a toroidal component; (ii) in the subsequent nonlinear reorganization of the system, ~ 90% 
of magnetic energy is lost in few Alfven timescales mainly through electromagnetic emission, and further de- 
creases on a much longer timescale; (iii) all stellar models tend to achieve a significant amount of magnetic 
helicity and the equipartition of energy between poloidal and toroidal magnetic fields, and evolve to a new 
configuration which does not show a subsequent instability on dynamical or Alfven timescales; (iv) the elec- 
tromagnetic emission matches the duration of the initial burst in luminosity observed in giant flares, giving 
support to the internal rearrangement scenario; (v) only a small fraction of the energy released during the pro- 
cess is converted into /-mode oscillations and in the consequent gravitational-wave emission, thus resulting in 
very low chances of detecting this signal with present and near future ground based detectors. 

Subject headings: stars: neutron 
merical 



gravitational waves — magnetohydrodynamics (MHD) — methods: nu- 



1. INTRODUCTION 

Neutron stars (NSs) are endowed with very intense, long- 
lived, large-scale magnetic fields, reaching strengths which 
are estimated to be of the order of 10 13 G at the magnetic 
pole for ordinary NSs, and around 10 15 G in the case of 
magnetars. Such extreme magnetic fields play a crucial role 
in the physics of NSs, affecting their structure and evolu- 
tion. They are involved in the processes through which NSs 
are observed, like the pulsar magnetic dipole radiation and 
the magnetically-powered burst activity of magnetars, and 
they have been recently recognised as essential to explain 
the quasi-periodic oscillations detected in the afte rmath of 
magnetar giant flares [see, e.g. . iGabler et al.l (120121) and ref- 
erences therein]. Moreover, they are responsible for defor- 
mations whic h may cause a significant emission of gravita- 
tional waves dBonazzola & Gourgoulhonl l 19961: ICutler 2002) 
and precession (Wasserman 20 Q3|). the y influence the thermal 
evolution of the star (IPons et al.ll2009l) . to list a few. 

All these processes depend on the magnetic field con- 
figuration inside the NSs, whose geometry is basically un- 
known. From observations of the spindown in pulsars the 
exterior magnetic field appears to be purely poloidal and 
mainly dipolar, but substantially different internal geometries 
can reproduce this external appearance. The importance of 
obtaining such information has motivated a significant ef- 
fort in studying possible equilibrium models of magnetized 
NSs, at first with simple geometries, e.g. purely poloidal 
or purely toroidal fields, and recently with mixed poloidal- 
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toroidal fields. The latest models, b uilt in Newtonian and 
genera l-relativistic framework, in clude l Tomimura & Eriguchl 
d200l: lLander & Jones! (120091 120121) : ICiolfi et alJ (l2009f 
120101) : iFuiisawa et al.l (120121) . where the so-called 'twisted- 
torus' configuration is considered. This particular geometry 
has been found as a result of the evolution of initial random 
fields in Newtonian magnetohydrodyn amic (MHD) simula- 
tions by lBraithwaite & Nordlundl (|2006h . 

Once a magnetic-field geometry is chosen, building a corre- 
sponding equilibrium configuration is not sufficient to assess 
whether this represents a good description of the NS interior. 
The magnetic field, in fact, should also be long-lived and thus 
stable on timescales which are much longer than the dynam- 
ical timescale. Assessing the stability of a given magnetic 
field configuration is not trivial and most of the work done 
on the subject concerns simple field geometries and nonro- 
tating stars. Until very recently, the problem has been only 
addressed with a perturbative analytic approach. These cal- 
culations established that both a purely poloidal field and a 
purely toroidal field are unstable in nonrotating stars, giving 
important predictions about the onset of the instability, but 
they could not predict the following evolution of the system. 
Only recently, thanks to the progress in numerical simula- 
tions, it has become possible to study these hydromagnetic 
instabilities by performing the fully three-dimensional (3D) 
MHD evolutions of magnetized relativistic stars. These simu- 
lations represent a very powerful tool, allowing to confirm the 
predicted features of the instability and to obtain information 
about the nonlinear dynamics of the process. In addition, the 
end-state of simulations can provide important hints about the 
preferred magnetic field configuration in magnetized stars. 

There is an additional and important motivation for study- 
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ing hydromagnetic instabilities in NSs. The induced global 
rearrangement of magnetic fields is a violent, strongly dy- 
namical p rocess, and soon after the m agnetar model was 
proposed dDuncan & Thompson! 1 1992b . this kind of pro- 
cess was suggested as a trigger m echanism of giant flares 
(iThompson & Duncanlll995L 12001b . Nowadays, this internal 
rearrangement scenario represents one of the leading mod- 
els to explain the phenomenology observed in magnetars, the 
other one involving a large-scale rearrangement of magnetic 
fields in the magnetosph ere surrounding the star dLyutikovl 
I2006t iGill & Heyj 12010b . Since both mechanisms can be 
present in a giant flare, the main question becomes whether 
most of the magnetic energy powering the flare is stored in- 
side the star or in its exterior magnetosphere. The basic 
tests on these models rely on the comparison of the predicted 
timescales and energies involved with the giant flare observa- 
tions. Hence, determining self-consistently the dynamics as- 
sociated to this kind of instability can provide important hints 
on the underlying mechanism. 

Moreover, magnetar flares (and in particular giant flares) 
are likely to be accompanied by a significant excitation of NS 
oscillations, in particular in the /-mode, which can then lead 
to a strong emission of gravitational waves (GWs). This pos- 
sibility has motivated recent searches for GWs in connection 
to magnetar flare s, published by the L IGO and Virgo collab- 
oration [see, e.g. lAbadie et al.l ( 1201 lb l. Semi-analytic efforts 
have been devoted to establishing the maximum amount of 
magnetic energy released in a magnetar flare, which, in turn, 
provides an upper limit on the energy emitted in GWs (llokal 
l200ltlCorsi & Owenll20Tlh . These upper limits are based on 
analytical calculations and simplified models, and can only 
provide rough, order-of-magnitude estimates. The assump- 
tion that all the available energy (which is at most of the 
order of the total magnetic energy) is converted into GWs, 
leads to the optimistic conclusion that the signal would be 
detectable with the n ext-generation ground-based detectors 
dCorsi & Owed 1201 lb. T his result has been questioned in 
ILevin & van Hoved d201 lb . where a simple perturbative anal- 
ysis is employed to show that only a small fraction of the mag- 
netic energy involved in a giant flare event can be actually be 
converted into /-mode oscillations and that the consequent 
GW emission is expected to be very weak. Again, referring 
to the internal magnetic field rearrangement scenario of giant 
flares, MHD simulations of hydromagnetic instabilities can 
provide a realistic picture of the GW signal produced, together 
with estimates of the fraction of energy which can be pumped 
into the /-mode and of the signal detectability. 

In this work we focus on the instability of purely poloidal 
fields, with the goal of shedding some light on all of 
the point s made above. Two parallel wor ks back in 
the '70s dMarkev & Tavleri [19731 IWrighj [l973h found that 
poloidal fields suffer from the so-called "Tayler" or "kink" 
instability, which manifests itself firstly in the neighbour- 
hood of the neutral line. This instability was recently studied 
with Newtonian numer ical simulations in the linear regime 
dLander & Jones! 1201 lb . or with non linear evolutions for a 
simpl ified model of newly born NS dGeppert & Rheinhardt 
2006b . and in the case of main-sequence stars dBraithwaite 
2007). The first 3D general-relativistic MHD simulations of 
the poloidal field instability i n NSs were presen t ed only last 
year, in two parallel works dLaskv et al.l 1201 it ICiolfi et all 
1201 lb . These studies reported similar results, despite some 
substantial difference in the approach (in particular in the evo- 
lution of magnetic fields outside the star), essentially confirm- 



ing all of the analytic predictions on the instability, and pro- 
viding some first hints about the nonlinea r rearrangement of 
the system. In addition. ICiolfi et ail (1201 ll) presented the first 
examples of gravitational waveforms triggered by the instabil- 
ity , which were subs equently considered more sy stematically 
by Zink et al.l (1201 lh for nonrotating stars and by lLasky et al.l 
d2012l) for r otating NSs. In th is paper we extend the work 
presented in ICiolfi et ail d201 lb . presenting additional infor- 
mation on the numerical infrastructure used and consider- 
ing the instability-driven evolution of a series of nonrotating 
NSs endowed with purely poloidal magnetic fields of different 
strength. 

The organization of the paper is as follows. In Sect. [2] we 
reconsider the setup of the system, improving in particular our 
treatment of the atmosphere. Within the new setup, we con- 
firm that our evolutions meet all the expectations on the onset 
of the instability , in a greement with the previous perturbative 
studies. In Sect. l3.ll we examine in more detail the nonlinear 
rearrangement of the system, performing much longer sim- 
ulations and gaining new subs tanti al insight on the final state 
reached by the system. Section [3~l2"l is dedicated to a discussion 
of the implications for the most-likely magnetic-field configu- 
rations in magnetized NSs and of the role played by magnetic 
helicity. We also study the emission properties of the system, 
relevant for the internal field rearrangement scenario of giant 
flares, estimating in Sect. I3.3l the timescale of the process and 
its electromagnetic luminosity , and discussing the detectabil- 
ity of the GW signal in Sect. 13.41 Both for electromagnetic 
and GW emissions, our estimates rely on a good agreement 
with the dependence on the magnetic field strength expected 
in the perturbative limit of weak magnetic fields, which al- 
lows us to extrapolate our results to lower and more realistic 
values than those actually considered in the simulations. Our 
conclusions are finally presented in Sect. [4] Unless specified 
differently, we adopt units in which c = 1, G = 1. 

2. PHYSICAL SYSTEM AND NUMERICAL SETUP 

Our physical system of interest is represented by a nonro- 
tating isolated neutron star surrounded by vacuum, initially 
endowed with a purely poloidal magnetic field permeating 
the star and extending to the exterior. The initial configura- 
tions are fully-relativistic self-consistent solutions generated 
with the multi-domain spectral-method code LORENE, de- 
velop ed at the Observatoire de Paris-Meudon dBocquet et alJ 
11995b . The stars are modeled as composed of a barotropic 
fluid obeying a polytropic equation of state p = Kp r , with 
T = 2 and K = 100. The reference unmagnetized star 
has a gravitational mass of 1.4 and a radius of 12.2 
km. Magnetic field strengths at the pole vary in the range 3 
Bp = Bx.o - £9.5 = (1 - 9.5) x 10 16 G. A stronger mag- 
netic field shortens the evolution timescale of the system, and 
the above choice makes our simulations computationally fea- 
sible. On the other hand, as we discuss in Sect. [3] most of our 
results can be extrapolated back to smaller (and more realis- 
tic) magnetic field strengths. 

We perform fully 3D general-relativistic MHD simulations 
of the system adopting the Cowling approximation, i.e. ne- 
glecting changes in the spacetime metric. Evolutions are ob- 
tained with the WhiskyMHD code. Most of the details on 
our mathematical and nu m erical setup are discussed i n dept h 
in iPollnev et all d2007b : iGiacomazzo & Rezzolla! d2007b : 

3 For simplicity, hereafter we will use a more compact notation for 
magnetic-field strengths in scales of 10 16 G, i.e. B N = N X 10 16 G. 
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iGiacomazzo et al.l (120091) : rGiacomazzo et al.l d201 ll) . One as- 
pect worth stressing is that, to guarantee the divergence-free 
character of the MHD equations, we use as an evolution vari- 
able the v ector potential inst e ad of the magnetic field, as de- 
scribed in IGiacomazzo et ail (1201 11) . Our standard grid setup 
consists of three refin ement levels using the Carpet driver 
dSchnetter et al.ll2004l) . with the finest one covering the entire 
star and having resolution h/M@ = 0.17 (~ 250 m). Our 
computational box extends to ± 54 M Q ~ 80 km and we im- 
pose no symmetries. 

To shorten the time for the development of the instability 
and thus reduce computational costs, a small perturbation is 
added to the initial velocity of the fluid. In particular, we add 
a ^-component of fluid velocity in the region surrounding the 
neutral line, where the instability is expected to be triggered, 
and with m = 2 azimuthal distribution. The strength of the 
perturbation corresponds to relative changes in the magnetic 
field of the order of 10~ 3 . We have verified that the instability 
occurs even in absence of initial perturbations and that no ap- 
preciable differences (besides that of reducing computational 
costs) are introduced in the dynamics by the perturbation. 

In our numerical setup, special care is paid to the treat- 
ment of the atmosphere. We recall that, as customary in 
finite-volume relativistic hydrodynamics simulations, our star 
is surrounded by a fluid at much lower densities (i.e. the "at- 
mosphere"), o btained by imposing a mi nimum rest-mass de n- 
sity [see, e.g. iBaiotti et all d2005h and iBaiotti et alj d2008l) l. 
In the atmospheric region, where the density is equal to the 
imposed floor value, we set the fluid velocity to zero to avoid 
the spurious accretion of the atmosphere onto the star. Within 
the common assumption of ideal MHD, also adopted in the 
WhiskyMHD code, this prescription would imply that mag- 
netic fields do not evolve in the atmosphere, which repre- 
sents a serious limitation from both the physical and numer- 
ical point of view. In particular, this rapidly leads to errors 
at the stellar surface and, for the magnetic field strengths we 
consider, ends the simulations prematurely. 

A self-consistent solution to this problem would require 
the implementation in general-relativity of the equations 
of resistive MHD, along the lines of the work reported 
by lPalenzuela et all d2009h and which has se en recent progress 
in the work of iDyonisopolou et al.l d2012h . Lacking for the 
time being this more systematic solution, a reasonable first 
approximation can be obtained by adding a resistive term to 
the evolution equation for the vector potential, i.e. 

d t A = v x B + i]AA, (1) 

where it = av — (3 fsee lAnton et al.l ( 120061) for the notation], 
77 is the scalar resistivity and, for simplicity, we take A to 
represent the Laplacian in a flat spacetime. We impose that 
resistivity is zero inside the star (thus retaining an ideal-MHD 
behaviour), up to a transition layer, where it continuously in- 
creases towards the atmospheric value 7/0 . This is shown in 
Fig. Q] where we plot the initial resistivity profile along the 
radial direction for our fiducial simulation with B p = i?6.5- 
The plot shows two different radial profiles, where we change 
the width of the transition region inside the star, whose sur- 
face r = R is indicated with a vertical solid line. As we will 
discuss below, the two choices have a different impact on the 
overall results. The resistivity profile is set to be an explicit 
function of the rest-mass density, namely f](p)/rio = g(p), 
where g is a Fermi-like function, with g = 1 in the atmo- 
sphere and g = inside the stellar core. In this way, the 
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Figure 1. Resistivity profiles along the radial direction for our fiducial simu- 
lation with Bp = -Bs.5- The two lines refer to the choice of a wide resistive 
transition layer (r ed dashed line) and corresponding to the one adopted in 
ICiolfi et al] 1201 lh . or of a thin layer (blue solid line). The vertical line marks 
the radius of the star R. 

resistivity will mimic the evolution of the rest-mass density, 
following the stellar surface in its evolution. 

Since the velocity is set to zero outside the star, eq. (fl]i re- 
duces there to dtA = 770 A A. A few remarks are worth doing 
about this limit and we start by considering the special rel- 
ativistic case for simplicity. Maxwell's equations in vacuum 
are simple wave equations, e.g. d^B = c 2 AB, so that a non- 
zero Laplacian of the magnetic field would be simply radi- 
ated away in electromagnetic waves. In the quasi-static limit, 
that is, in the limit in which the timescale associated with 
magnetic-field variations is large compared to the light travel 
time L/c (where L is the typical lenghtscale of field varia- 
tions), any episodic time variation in the magnetic field pro- 
duced by the dynamics in the stellar interior would be rapidly 
radiated away, leading to an exterior magnetic field which is 
again with AB ~ 0. In our system, where the quasi-static 
limit represents a good approximation, a similar behaviour 
outside the star can be obtained by evolving the magnetic field 
according to a diffusion equation 4 , dtB = 770 AB. In this 
case, the non-zero Laplacian components of B are removed 
not through the propagation of electromagnetic waves, but 
through resistive dissipation. Hence, as long as one is not in- 
terested in the precise dynamics of the exterior magnetic field, 
the two recipes are effectively equivalent and the addition of 
a resistive term in eq. ([]]) has the effect of mimicking the be- 
haviour of the Maxwell equations in the vacuum outside the 
star, at least within the quasi-static approximation. Indeed, 
this approach is not novel, but it has already been consid- 
ered in the literature, for instance, by lBraithwaite & Nordlundl 
(2006). 

This logic is valid as long as the external evolution takes 
place in the quasi-static limit, and we fulfil this requirement 
by suitably choosing the timescale of magnetic diffusion via 
the choice of the resistivity 770. More precisely, we set ?yo to 
be always high enough (hence with sufficiently short associ- 
ated resistive timescales) so as to keep the exterior of the star 
always close to the condition AB ~ 0. Although only an ap- 
proximation, the approach discussed above allows us to have 
a dynamical exterior magnetic field, which adjusts itself to the 
changes triggered in the interior by the development of the in- 
stability. The price to pay is that, in addition to the one lost 
because of the numerical resistivity, the magnetic energy of 

4 Note that the same equation holds for the vector potential, dtA = 
770 A A 
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Figure 2. Poloidal field instability in our fiducial simulation with B p = -E?6.5- Panels show different stages of the evolution (from left to right t = 3.5, 7.5, 
60 ms) with meridional view (top row) and equatorial view (bottom row). Shown with vector lines are the (global) magnetic-field lines, while the colors show 
the intensity (in Gauss) of the toroidal magnetic field only; also reported are three rest-mass isodensity contours near the stellar surface, corresponding to 
p = (0.02, 0.2, 2) X 10 13 g/cm 3 . 



the system is not conserved. Per-se, this would not be particu- 
larly problematic, since in a realistic system one would expect 
the magnetic energy to be radiated away, but it does introduce 
a dependence of the problem on the profile chosen for the re- 
sistivity, since different profiles will be responsible for energy 
losses. 

As we will discuss in more detail the following Section, 
in fact, during the initial development of the instability, the 
magnetic -field modifications are still very small compared to 
the background field, and the system is not expected to show 
visible changes, with the magnetic energy being essentially 
conserved. Yet, because of the resistive layer inside the star 
introduced via eq. (fl]i, magnetic-energy losses will be present 
even in the early stages of the evolution, simply because the 
magnetic field will be divergence free but not Laplacian-free 
in the outer layer of the star. 

In our previous work dCiolfi et al.H2"oTll) . we had chosen a 
rather wide resistive transition layer (cf. red dashed line in 
Fig. [TJ, which resulted in significant magnetic-energy losses 
from the very beginning of the simulation. To minimize 
these losses, we have here considered a thinner resistive tran- 
sition layer (cf. blue solid line in Fig. [TJ, thus involving a 
smaller stellar volume of the star where resistivity can act. As 
a result, the magnetic-energy losses during the exponential- 
growth phase (i.e. during the first ~ 3.5 ms) go from ~ 30% 
in the case of a wide layer, to below ~ 2% in the case of 
a thin layer. In the subsequent stages of the instability, the 
differences between the two prescriptions are much smaller 
and this is because the dynamics of the field is much less 
influenced by the properties of the resistive layer at the sur- 
face of the star. Yet, since we have shown that different lay- 
ers lead to different energy losses, it is reasonable to wonder 
whether the properties of the transition region can be impor- 



tant also for the subsequent evolution of the system. As we 
will argue in detail in the Appendix, results are effectively 
independent of 770 as long as a suitable value for i]o is cho- 
sen for the different magnetic field strengths; more specifi- 
cally, if the resistivity is chosen to scale linearly with B p as 
77o/M Q = rfo x (B p /B 6 . 5 ) = 0.1 x {B p /B 6 . 5 ) (see Ap- 
pendix). 

3. RESULTS 

3.1. Poloidal field instability 

We next describe the evolution of the system for our fiducial 
set of initial models, from the onset of the instability to the 
nonlinear rearrangement of magnetic fields. First of all, we 
briefly recall the basic expectations of the perturbative stud- 
ies on the onset of the instabi lity for a purely poloidal field in 
nonrptating magnetized stars (Ma rkev & Tavlerill973fclw"rightj 
1 1973b : (i) the instability first develops in the region of closed- 
field lines surrounding the neutral line; (ii) toroidal magnetic 
fields are generated in this region and grow exponentially until 
their local intensity is comparable to the poloidal one, which 
corresponds to the saturation of the instability; (iii) the insta- 
bility saturates in about one Alfven time 5 ; (iv) the timescale 
associated with the exponential growth of the toroidal field 
scales as 1/B with the magnetic field strength. As we will 
show in the following, our simulations meet all of these ex- 
pectations. 

A first visual overview of the dynamics of the system is 
given in Fig. [2] where we present snapshots of the star and its 
magnetic field in the meridional (y, z) (top row) and equato- 
rial (x, y) planes (bottom row) at three different stages of the 

5 An estimate of the Alfven time is given by ta ~ ?>R%J 47r(p) /-Bp, 
where (p) is the average rest-mass density. Hence, ta ~ 3 ms for our 
fiducial model with B p = Bq. 5. 
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Figure 3. Left Panel: Evolution of total magnetic energy normalized to the initial value, for different initial magnetic field strengths: i?6.5 (continuous black 
line), -B5.0 (dashed red line), -B3.5 (long-dashed blue line) and -B1.5 (dot-dashed magenta line). Middle panel: Evolution of poloidal and toroidal magnetic 
energies normalized to the initial total magnetic energy (log scale), for Be. 5 , -Bl.s. Right panel: Ratio of toroidal and poloidal magnetic energies versus time, 
for the same collection of simulations shown in the left panel. 



evolution. In this example B p = £?6.5- The left panels refer 
to the time of saturation of the instability, when the magnetic 
field starts to show visible modifications (t = 3.5 ms). From 
the equatorial view we note that the initial axisymmetric ge- 
ometry is lost and replaced by a non-axisymmetric structure. 
Around 7.5 ms (central panels) the instability has fully de- 
veloped. As expected, the closed line region is filled with 
toroidal fields of strength comparable to the background, re- 
sulting in vortex-like structures of magnetic field lines. At this 
point the initial geometry of the system will be gradually lost. 
This stage is the most dynamical and the rapid modifications 
of the field lead to the expulsion of matter from the surface of 
the star (~2x A/© in this example). The last panels re- 
fer to the end state of our simulation (t = 60 ms). At this stage 
there is no trace of the initial geometry and the system has lost 
most of its magnetic energy (see discu ssion below ) . In c om- 
parison with the snapshots presented in lCiolfi et al.l (1201 ll) for 
the same magnetic field strength, the evolution appears more 
violent. The significant magnetic energy loss due to resistiv- 
ity in a thicker transition layer, in fact, had the effect of re- 
straining the dynamics, resulting in a less dramatic and more 
ordered evolution. 

In Fig. [3] we plot the evolution of the poloidal, toroidal and 
total magnetic energies. In practice, the evolution of the total 
magnetic energy (left panel) provides information on the var- 
ious stages of the evolution: the initial phase up accompanied 
by very small variations, the sudden drop of magnetic energy 
with a loss of ~ 90% of magnetic energy in few ms to tens 
of ms depending on the initial magnetic field strength, and the 
final stage with a slower evolution as the energy decreases to 
few percent of the initial value and the system reaches a quasi- 
stationary equilibrium 6 . Assessing the equilibrium properties 
of the new configuration is complicated by the fact that, by 
construction, our system would suffer from resistive losses 
due to a thin layer inside the star and outside of it, even if a 
stable equilibrium had been reached. Hence, these residual 
losses make it difficult to determine unambiguously whether 
the new, post-instability configuration is a stationary one or 
not. What is evident is that the system's properties are not 

6 Note that we refer to a "quasi-stationary equilibrium" and not to a stable- 
equilibrium because, strictly speaking, the latter is impossible to prove with 
numerical simulations over a finite amount of time. 



changing significantly (e.g. the rest-mass density) or, if chang- 
ing, they are doing so much more slowly than during the in- 
stability (e.g. the magnetic field). We interpret this behaviour 
as (partial) evidence that the new magnetic-field configuration 
has reached a quasi-stationary state or that, if still intrinsically 
unstable, the growth time of the instability is much larger than 
the one we can possibly investigate. As we will comment be- 
low, this conclusion is further supported by the dynamics of 
the rest-mass density. 

The overall evolution does not change with smaller initial 
field strengths, apart from the expected increase of timescales. 
Figure [3] also reports the poloidal and toroidal field ener- 
gies in a logarithmic scale (middle panel). Note that the the 
toroidal energy grows exponentially in the initial stages of 
the evolution and that the growth rate depends on the ini- 
tial magnetic-field strength, the slope being steeper for larger 
magnetic fields. The poloidal component, on the other hand, 
remains almost unchanged during this initial stage, until the 
nonlinear rearrangement of the magnetic field starts. Then, 
the system evolves losing most of the poloidal energy, while 
the toroidal energy experiences smaller variations, thus result- 
ing in a growing ratio of toroidal and poloidal energies. In the 
last phase of the evolution, the ratio of the poloidal-to-toroidal 
magnetic energies stabilizes in the range ~ 0.7 — 1, indicat- 
ing that at this stage the toroidal magnetic field provides a dy- 
namically important contribution to the final quasi-stationary 
balance (right panel of Fig . [3} . 

From a closer look at the details of the dynamics, we notice 
a qualitative difference between stronger and weaker mag- 
netic fields, the latter having a much smoother evolution. For 
Bp = B\x) — B2.5 the system undergoes indeed less dra- 
matic modifications, where the initial magnetic field geome- 
try is partially maintained. This is evident in Fig. |4] where 
we show snapshots of the evolution of a star with initial mag- 
netic field strength B p = B1.5. These differences are relevant 
when we consider the emission properties of the system and, 
in particular, the emission of gravitational waves. 

In Fig.|5]we show instead the growth time r of the toroidal 
field energy in the exponential phase of the evolution, as a 
function of the initial magnetic field strength. The predicted 
linear scaling is very well satisfied in the full range of mag- 
netic field strengths considered (the red dashed line represents 
a linear fit to the data). 
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3.2. Final magnetic field configuration 

Since the new simulations have been carried out on 
ti mescales which are much longer than the ones investigated 
in ICiolfi et all (IMTTh . we are in a much better position to 
discuss the properties of the final magnetic-field configura- 
tion. Also important in this context is the evolution of purely- 
hydrodynamical quantities, which can provide important ad- 
ditional information on the development of the instability 
and on the quasi-stationary state approached by the system 
in the final stages. One of such quantities is the central 
(i.e. maximum) rest-mass density. We recall that the initial 
purely poloidal magnetic field is responsible for a quadrupo- 
lar deformation of the star, whose resulting shape is slightly 
oblate (i.e. with a positive ellipticity e), despite being non- 
rotating. This deformation lowers the central rest-mass den- 
sity of the star with respect to the unmagnetized case, and 
the difference is of the order of Ap max / p max (B = 0) = 
1 — Pmax 1 1 Pmax{B = 0) ~ e. Since the ellipticity scales as 
e ~ B 2 , when the system loses a factor ~ 10 2 in poloidal 
magnetic energy, Ap max should be also reduced by about the 
same factor. In addition, if there is a significant toroidal com- 
ponent, this tends to deform the star in a prolate shape, con- 
trasting the effect of the poloidal field. As a consequence, in 
our simulations we expect the final Ap max to be much smaller 
than the initial one, i.e. the star basically recovers the spheri- 
cal shape of the unmagnetized case. 

In Fig.|6]we show the evolution of the central rest-mass den- 
sity for different initial magnetic field strengths. Note that our 
different models are all built with the same initial central den- 
sity, which gives a percent difference in the central density of 
the unmagnetized star corresponding to each model. Clearly, 
also for the rest-mass density the evolution is much more dra- 
matic for stronger initial magnetic fields, with excursions in 
the central density that become larger and more rapid as the 



initial magnetic field strength is increased. In all cases, how- 
ever, the central density mimics the evolution of the magnetic 
energy and reaches an approximately constant value, with an 
overall jump which scales roughly as ~ B 2 . The absence of 
additional evolution in the rest-mass density after the instabil- 
ity has fully developed provides the indication that the new 
configuration is, at least hydrodynamically, stable, and that, if 
the new magnetic-field configuration is about to lead to a new 
unstable evolution, it will do so on much larger timescales. 
For initial magnetic field strengths of 6.5 and 5 x 10 16 G, the 
change is so violent and rapid that it leads to high-frequency 
oscillations in the central density, while we do not observe the 
same effe ct with smaller fields. As we will further discuss in 
Sect. 13.41 a rapid variation of central density can also affect 
the emission properties of the star, e.g. introducing a signifi- 
cant modulation of the gravitational-wave signal. 

Other useful information on the final state can be obtained 
by looking at the evolution of the total magnetic helicity H m . 
This quantity is associated with the topological properties 
of a given magnetic field geometry, a nd measures t he de- 
gree of linkage of magnetic field lines (Moffatt 1969). The 
initial purely poloidal configuration has H m = because 
only regions with a mixed poloidal-toroidal field contribute 
to the magnetic helicity 7 . We already noticed that the ratio of 
toroidal and poloidal energies grows in time and then tends to 
stabilize to values in the range ~ 0.7 — 1. The total magnetic 
helicity, even if related to the energy in the two magnetic-field 
components, encodes an independent information about the 
topology of the field, and it is therefore useful to investigate 
its evolution in our system. 

We recall that although in ideal MHD the total magnetic he- 

7 We note, however, that even a mixed field can be constructed such that 
H m = 0, e.g. in axisymmetry, through the superposition of a poloidal mag- 
netic field and of two toroidal ones of opposite polarity. 
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Figure 5. Instability growth time as a function of the initial magnetic field 
strength (black circles). The red dashed line represents a fitted linear scaling 
as expected from perturbation theory. 

licity is conserved ( IWoltierll 19581) . a highly-conducting fluid 
star in vacuum with a thin resistive layer does not represent 
an ideal-MHD system, and magnetic helicity conservation is 
therefore not guaranteed [of course helicity is also produced 
because of the intrinsic nonzero numerical resistivity, but the 
latter is much smaller than the one introduced in eq. Q)]. We 
compute the total magnetic helicity as 

H m = f H^y/=g<fx, (2) 

where = *F a P Ap is the magnetic helicity 4-current 
and Tj t is the spa tial hypersurface at a given time t [see 
iBekensteTnl (119871) and references therein]. Since our ini- 
tial helicity is H m (t = 0) = 0, we normalize H m to the 
helicity of an axisymmetric twisted-torus configuration in 
which the poloidal field is arranged as in our initial condition 
and the toroidal field, confined in the region of closed field 
lines around the neutral line, is uniform and with an energy 
equal to the poloidal one, i.e. E po i = E tor = E m /2. Indicat- 
ing with H m such reference helicity, we can estimate it as 
\H m \ ~ r N ^Epoi E tor = r N E m /2, where r N (-8 km) is 
the radius of the neutral line and E m is set to coincide with the 
magnetic energy of the stellar model under exam. Note that 
we are not concerned here with the sign of H m (or of H m ), 
since the latter depends on whether toroidal fields are directed 
along the positive or negative (^-direction and a global trans- 
formation — ^ — would not change the properties of the 
system. 

We can now monitor the evolution of H m . As long as 
\H m /H m \ -C 1 the total magnetic helicity is small, while 
\H m /H m \ ~ 1 is an indication that there is a significant 
magnetic helicity. In Fig. [7] (top panel) we show the evolu- 
tion of the normalized helicity for the different initial mag- 
netic field strengths. Note that in all cases \H m /H m \ grows 
in time, and around 60 ms it has reached a substantial frac- 
tion of 1. This fraction is higher for weaker fields, where 
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Figure 6. Evolution of central rest-mass density for different initial magnetic 
field strengths. From top: 6.5 (black), 5 (red), 3.5 (blue), 1.5 X 10 16 G 
(magenta). 

the evolution is smoother, i.e. \H m /H m \ ~ 0.6 for -B1.5, 
while \H m /H m \ ^0.3 for B&.§. This behaviour supports the 
idea that configurations with a significant amount of magnetic 
helicity are more stable, as a l ready suggested in the litera- 
ture rsee lBraithwaite & Spruitl (120061) and references therein]. 
Such indication is compatible with the observed instability of 
the initial purely poloidal field and the production of a toroidal 
component, which could be interpreted as an attempt of the 
system to produce magnetic helicity. 

We have seen that the system can achieve a significant 
amount of magnetic helicity, but in the meanwhile most of 
its magnetic energy is lost. In order to weigh the absolute 
variation of magnetic helicity, we can compare H m with the 
initial value of H m (bottom panel of Fig . O . The ratio in this 
case is much smaller, < 1 — 3%, indicating that the magnetic 
helicity produced would represent a small amount for the ini- 
tial star, and becomes significant only because the system has 
lost most of its magnetic energy. 

On the basis of these results and bearing in mind the lim- 
its of our approach, we conjecture that magnetic helicity could 
play an important role as a stabilizing element for a fluid mag- 
netized star. If this is the case, we can give a natural interpre- 
tation of the evolution of our system, (i) A purely poloidal 
magnetic field is unstable and a significant modification of its 
topology is necessary for stabilization, (ii) While the system 
is evolving the rate of change of magnetic helicity is small 
compared to the one of magnetic energy, thus the system has 
to lose most of its magnetic energy in order to significantly 
alter its topology and reach a more stable state. 

As a complement to the information obtained from our non- 
linear relativistic calculations, we should remark that the pos- 
sibility of a stable equilib rium has been questioned in the case 
of a b arotropic fluid star dReisenegge r 200l lLander& Jones! 
2012). A barotropic equation of state does not take into 
account the effect of stable stratification associated to com- 
position gradients, which could play an important role in 
determining the stability of magnetic equilibria in relativis- 
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Figure 7. Top: Evolution of magnetic helicity normalized to H m (see text), 
for different initial magnetic field strengths: 6.5 (continuous black line), 3.5 
(long-dashed blue line) and 1.5 X 10 16 G (dot-dashed magenta line). Bottom: 
Evolution of magnetic helicity normalized to the initial value of H m , for the 
same set of simulations. 

tic NSs dReisenegger & Goldreichlll992t lReiseneggerj|2009l) . 
Stable stratification, indeed, offers an additional degree of 
freedom which favours the hydrostatic balance of fluid and 
magnetic forces. A strong indication in support of stable 
equilibria in stratified magnetized stars has been obtained for 
main-sequence s tars, where stratification is pro vided by en- 
tropy gradients dBraithwaite & Nordlundl 120061) . while there 
is no such evidence for barotropic NSs. 

3.3. Electromagnetic emission 

As discussed in the Introduction, one of the leading expla- 
nations of magnetar giant flares assumes the event as due to 
a sudden, global rearrangement of the internal magnetic field, 
i.e. the same kind of process that we are studying. A giant 
flare starts with an initial burst of 0.1 — 0.5 s duration, fol- 
lowed by a long pulsating tail lasting hundreds of seconds. 
In such scenario, the energy emitted in the initial burst is a 
significant fraction of the total magnetic energy of the star, 
while the duration of the burst is dictated by the timescale 
of the internal field rearrangement. We have also remarked 
repeatedly that in all cases considered, our magnetized stars 
lose most of the initial magnetic energy within few tens of ms. 
This energy is converted in very small part into fluid motions, 
e.g. /-mode oscillations of the star with consequent emission 
of gravitational waves (see Sect. 13.4b . while most of it is dis- 
sipated in the resistive layer inside the star and then in the 
atmosphere, thus mimicking the radiative losses that would 
be measured if the stars were actually in vacuum. We can 
therefore use the information about the dissipated magnetic 
energy to deduce, within the approximation of our approach, 
an order-of-magnitude estimate of the electromagnetic lumi- 
nosity L em associated with the hydromagnetic instability. In 
addition, we can compare the luminosity and duration of the 
emission computed in this way with the observations of giant 
flares. 

In practice, most of the emission comes from the initial fast 



drop of magnetic energy by ~ 90%, corresponding to an ini- 
tial spike in L em . Taking that as a reference, we measure the 
duration of the spike and compute its peak and (time) average 
luminosities for the different initial magnetic field strengths 
considered. Note that the expected scaling for the luminos- 
ity is L em oc B^, since the timescale of the emission goes 
as l/B, while the magnetic energy scales as E m oc B 2 . Our 
estimate for the average luminosity, obtained by fitting for the 
different magnetic field strengths, gives 

(L em ) = 1.9 x (-^Tq) x 10 48 erg/s. (3) 

Similarly, for the reference value of a magnetar-type mag- 
netic field of Bp = 10 15 G, we obtain a duration of the ini- 
tial spike of ~ 0.7 s and a peak luminosity of L em pea k ~ 
5.3 x 10 48 erg/s. 

We can now compare with the observations. In the bright- 
est giant flare detected, the one in 2004 from SGR 1806- 
20, the initial spike lasted ~ 0.5 s, had a peak luminosity of 
~ (2 — 5) x 10 47 x (d/15 kpc) 2 erg/s, where d is the distance 
of the source, and accounting for the total isotropic energy ra- 
diated (which was ~ 99% of the total energy emitted in the 
whole giant flare event, including the ~ 400 s tail), the aver- 
age hirmnosity was~ (0.3 — 1) x 10 47 x (d/15 kpc) 2 erg/s 
[see lMereghettil ( 120081) and references therein]. The duration 
of the initial burst from SGR 1806-20 is clearly compara- 
ble with our estimate. The difference in average luminosity 
(our estimate is ~ 20 — 60 larger), then, essentially reflects 
a difference in the electromagnetic energy released, which 
does not constitute a significant limitation for the scenario: 
while in our system the magnetic energy is almost completely 
lost, it is perfectly reasonable that in a more realistic situation 
the magnetic field would migrate from one configuration to 
another with a smaller jump in magnetic energy. Note that 
the timescale of the process, on the other hand, is essentially 
controlled by the initial internal magnetic field strength and 
poorly depends on the overall jump in magnetic energy. Inter- 
estingly, also the ratio between the estimated peak and aver- 
age luminosities is in good agreement with the observations. 

If the comparison with the phenomenology of from SGR 
1806-20 leads to a reasonably good match in the duration and 
energetics, thus providing a substantial support to the inter- 
nal field rearrangement scenario, the comparison with other 
observations is not as striking. In particular, the average lu- 
minosity of the initial spike in the case of the other two gi- 
ant flares observed (1979 from SGR 0526-66 and 1998 SGR 
1900+14) was < 10 45 erg/s, with a comparable duration of 
~ 0.25,0.35 s. In these cases, therefore, the mismatch in 
the energy losses is far larger, although still acceptable if the 
field rearrangement is to a new configuration with comparable 
magnetic field strengths. 

In summary, although the dynamics of our system is over- 
simplified when compared to the complex phenomenology 
shown by a realistic magnetar giant flare, the overall agree- 
ment in the duration of the burst and in its energetics when 
compared with the observations of the giant flare from SGR 
1806-20, lends support to the suggestion that the basic phe- 
nomenology is that of an internal-field readjustment. To the 
best of our knowledge, this is the first time that a compar- 
ison between general-relativistic MHD simulations and the 
magnetar phenomenology has been attempted. Clearly, a lot 
more work needs to be done to improve our self-consistent but 
crude modelling. 
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Figure 8. Gravitational wave signal for different initial magnetic field strength B p , in the range B p = 
in X (+) polarization; d is the source distance, which needs to be specified to obtain the amplitudes hy 



— i?6.5- Top (bottom) panels refer to the amplitude 
Note that the y-scale changes with B p . 



3.4. Gravitational wave emission 

In the Introduction we have pointed out the potential rel- 
evance of our study for GW observations in connection to 
magnetar g i ant fla res. The basic idea, already presented in 
ICiolfi et all d201 lb . is that the instability triggers stellar oscil- 
lations at the /-mode frequency, with the consequent emission 
of GWs, similarly to what should happen in association with a 
magnetar giant flare. In what follows we provide a systematic 
assessment of this emission for the different initial magnetic 
field strengths considered and an estimate of the detectabil- 
ity of the GW signal for the planned advanced GW detectors, 
i.e. advanced LIGO and advanced Virgo, and the new g enera- 
tion ones, e.g. Einstein Telescope dPunturo et al.ll2010l) . 

The GW signals produced in our simulations are shown 
in the different panels of Fig. [HJ where we report the am- 
plitudes in the x and + polarizations for the various initial 
magnetic field strengths B p , ranging from £?i.q to i?6.5- Be- 
sides the differences in the overall amplitude, that naturally 
grows with B p (see also discussion below), the waveforms 
show also other differences with the magnetic field strength. 
These include variations in the transient stage of the instabil- 
ity and low-frequency modulations. For example, for i3g.5 we 
can observe a significant modulation at ~ 25 Hz, which could 
be associated with the jump in central density (cf. Fig. [6), 
whose timescale is compatible with the period of the mod- 
ulation. A similar modulation appears also for i?2. 5- Rather 
irregular transients characterize instead the initial amplitudes 



for B3.5 and B$x>. In both cases, in fact, we can notice sudden 
jumps in the signal, which correspond to the highly dynami- 
cal phases in the evolution of the system. This is indeed con- 
firmed by the evolution of the total magnetic energy in Fig.[5J 
where bumps appear around ~ 12 ms and ~ 22 ms for B3 .5 
and B5.0, respectively, while they are not present in the other 
cases. 

Despite these differences, all signals are essentially domi- 
nated by oscillations at the /-mode frequency, as it is clear 
from Fig. [9] where we report the Fourier transform of the GW 
amplitude in both the x and + polarization for B1.5 , B3.5 and 
-Bg.5- The /-mode we can read off the plot is about 1.85 kHz 
for the lowest field strength, with a shift to lower frequencies 
for higher magnetic fields (this shift is of ~ 8% for Bq^). We 
recall that since our evolutions are in the Cowling approxima- 
tion, the corresponding /-mode frequencies are notoriously 
larger of ~ 15 % than the correct ones computed in full gen- 
eral r elativity dDimmelmeier et aT]|2006t iTakami & Rezzollal 
201 1). This property will be taken into account when consid- 
ering the detectability of the signal. 

Given the complexity of the signal, with large rapid and sec- 
ular variations, the determination of the signal-to-noise ratio 
(SNR) cannot be done by simply looking at the maximum- 
minimum amplitudes, but rather by computing a strain which 
represents a suitable time average. This is indeed what can be 
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Figure 9. Fourier transform of h x and for B p = -B1.5 (lower magenta 
line), -B3 5 (intermediate blue line), and Be. 5 (upper black line), assuming a 
source distance of 10 kpc. The vertical line marks a frequency of 1850 Hz. 

done through the root-sum-square amplitude h laa , defined as 

1 1/2 



h(t) 2 dt 



(4) 



where h 2 (t) = h 2 x (t) + h\(t). The amplitudes in the 
two polarizations are computed considering a source within 
the Galaxy, i.e. at a fiducial distance of 10 kpc. If, as in 
our case, the signal is dominated by a single frequency /o, 
then the SNR can be estimated using the simple expression, 

S/N = h Ias / y/S h {fo), where y/S h (fo) is the strain-noise 
amplitude of the detector at the frequency f . When comput- 
ing h rsa we also need to specify the duration of the signal (in 
addition to the distance of the source). Typical estimates of 
the /-mode damping time t>, w r a nge between 100 m s and 
1 s (lAndersson & Kokkotasl 119981; iBenhar et al.l 120041) . and 
here we assume r GW = 100 ms. 

The resulting h rss amplitudes are reported in Fig. [TO] 
where we also show with dotted lines the strain-noise am- 
plitude of advanced LIGO and Virgo and of the Einstein 
Telescope, and where we have assumed fo ~ 1550 Hz 
to correct for the Cowling approximation. Furthermore, 
the "noise" is estimated to be \fShJjo) ~ 7.4 x 
10 _24 Hz -1 ^ 2 for advanced LIGO and and advanced Virgo, 
and yS h (fo) ~ l(r 24 H z~ 1/2 for the Einstein Telescope 
fsee lAndersson et aTl (1201 ll) . where the same strain-noise am- 
plitudes were considered] . 

Because in our simulations we have considered magnetic- 
field strengths that are at least an order of magnitude larger 
than those typically associated with magnetars, we need to 
rescale our results for h ISS to smaller field strengths in order to 
assess the GW signal that could be expected if the instability 
is associated to a magnetar giant flare. This is rather straight- 
forward to do since the expectat i on fro m perturbative analyses 
[see, e.g. , iLevin & van Hover] ( 1201 1|) 1 is that, independently 
on the mechanism considered, the amplitude of the excited f- 
mode by the flare should scale as the magnetic energy, i.e. as 



Figure 10. h VS s [y/s\ versus B p in Gauss assuming a distance of 10 kpc 
and a damping time of t gw = 100 ms. Horizontal lines mark the strain- 
noise amplitude of the considered GW detectors at the /-mode frequency fo. 
Finally, the red line is obtained by imposing a quadratic scaling (ex B 2 ) and 
fitting the results for the lowest magnetic fields (first four points). 

oc B 2 . This prediction is valid as long as the magnetic field 
strengths are sufficiently small so that the corresponding mag- 
netic energy can be considered as a small perturbation to the 
total binding energy. The results shown in Fig. [TO] confirm 
that a quadratic scaling is a rather good approximation to the 
computed data for B p < i?3.o- This is highlighted in FigfTUl 
by the red solid line, which represents a quadratic fit to the 
GW amplitude in the case of low-magnetic fields. We find it 
very reassuring that our data does show the expected scaling 
behaviour in the weak-field limit, which instead app ears to be 
absent in the scaling reported b y IZink et all d201 ll) and sub- 
sequently bv lLaskv et alj (120121) for rotating stars. Of course, 
this scaling is then broken as the magnetic field is increased 
and a steeper dependence is then found, which appears to sat- 
urate for Bp > .B5.0 G. 

This quadratic fit allows us therefore to extrapolate back 
to even smaller values of the magnetic fields and obtain the 
following predictions for the SNRs for the different detectors: 



S_ 

N 



AdvLIGO-Virgo 
S_ 

N ET 



1.2 x 10~ 4 x 



0.9 x 10" 3 x 




, (5) 



(6) 



We conclude that the GW signal produced by the /-mode os- 
cillations triggered by the hydromagnetic instability in asso- 
ciation with a magnetar giant flare will be undetectable for 
realistic values of the magnetic field, i.e. B p < 10 15 G and 
marginally detectable by third generation detectors for unre- 
alistic magnetic fields, B p > 2 x 10 16 G. Because a longer 
damping time, say of t ow = Is, would yield a gain of only a 
factor ylO, we do not expect that the prospects of a detection 
would improve if the /-mode oscillations would last a factor 
ten more in time. 
Understanding why the detectability of this GW signal is 
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so hard in practice can be made easier if we take into account 
how much energy is actually lost to gravitational waves and 
compare it with the amount of dissipated magnetic energy. 
More specifically, if h ISS is the root-sum-square GW ampli- 
tud e, the energy emitted in G Ws can be computed as [see, 
e.g. ILevin & van Hovea (1201 ll) l 

Egw = — Q — d h iaa , 

"»"" , i«(A)'«&)' t!l 

where we retained factors of c and G, and d is the source 
distance 8 . Comparing expression (|7]i with the corresponding 
estimates of the dissipated magnetic energy (O, it is easy to 
realize that the two energies differ by more than ten orders of 
m agnitude. Hence, assu ming Egw ~ E m as done recently 
bv lCorsi &Owenl d201 lh would indeed lead to a S/N ~ 30 
for advanced LIGO, but it would also be overly optimistic 
and rather unjustified on the basi s of our calculations, which 
instead confirm the expectation of ILevin & van Hovenl(l20TTl) 
that the energy actually converted into GWs is only a small 
fraction of the total energy available. As a result of this in- 
efficient conversion, the GW signal is too weak to be de- 
tected with the near future detectors. A similar conclusion 
was al so reached by ILevin & van Hoveri (1201 ll) : IZink et all 
d20llh . and more recently by lLaskv et al.l ( 120121) for rotating 
stars. 

4. SUMMARY AND CONCLUSION 

We have performed 3D general-relativistic MHD simula- 
tions of nonrotating magnetized NSs endowed with a purely 
poloidal magnetic field and studied the development of the hy- 
dromagnetic instability that develops dynamically. This work 
represents an exte nsion of our previous study on the subject 
(ICiolfi et al.l 1201 ll) . where we have improved our treatment 
of the atmosphere, drastically reducing the undesired energy 
losses due to the transition between the ideally conducting 
stellar interior and the resistive exterior, and we have per- 
formed simulations on much longer timescales, which have 
allowed us to gain essential information about the final con- 
figuration reached by the system. The resulting overall picture 
is much clearer and conclusive. 

As expected from perturbation analyses, the instability is 
first triggered in the region of close field lines and is accompa- 
nied by the production of a toroidal magnetic field. When the 
growth of the latter saturates in about one Alfven timescale, 
the toroidal field has reached a local strength which is com- 
parable to the poloidal one and the axisymmetry of the initial 
configuration is lost. At this point, the most dynamical phase 
of the nonlinear evolution takes place, with major modifica- 
tions of the magnetic field leading to a strong electromagnetic 
emission carrying away ~ 90% of the magnetic energy in few 
Alfven timescales. At the same time, a small fraction of mag- 
netic energy is converted into stellar oscillations, mostly at 
the /-mode frequency, which cause the emission of gravita- 
tional waves. The subsequent evolution proceeds on longer 
timescales, with further loss of magnetic energy. 

Since in the post-instability phase the magnetic field is con- 
tinuously changing, losing strength because of resistive dissi- 
pation, it is difficult to determine whether the corresponding 

8 Note that d 2 h^ BB does not depend on the source distance, while h TBB 
alone does. 



configuration is a stationary one or not. The only robust ev- 
idence is that the variations in the hydrodynamical and elec- 
tromagnetic quantities are much smaller than those during the 
instability, taking place on larger timescales. We therefore in- 
terpret this behaviour as evidence that the new magnetic-field 
configuration has reached a quasi-stationary state or that, if 
still intrinsically unstable, the growth time of the instability 
is much larger than the one that can be possibly investigated 
numerically. 

Because in this quasi-stationary configuration the ratio of 
the toroidal and poloidal magnetic energies tends to unity, 
and because the growth of toroidal magnetic field at the ex- 
pense of the poloidal one is also accompanied by an increase 
in the magnetic helicity, we are led to conclude that an equi- 
librium magnetic field configuration with a significant amount 
of magnetic helicity and comparable poloidal and toroidal 
magnetic field energies could be a preferred one for stabil- 
ity. Bearing this in mind, we remark that it is still unclear if 
stable equilibria exist for a simple barotropic fluid star, and 
that other stabilizing contributions, such as a stable stratifi- 
cation, may be necessary to obtain long-lived magnetic field 
configurations. 

The violent reorganization of magnetic fields induced by 
the development of a hydromagnetic instability has been pro- 
posed as a possible mechanism to explain giant flares in mag- 
netars. Using our simulations and in particular the informa- 
tion about the dissipated magnetic energy, we have deduced, 
within the approximation of our resistive approach, an esti- 
mate of the electromagnetic luminosity L em associated with 
the hydromagnetic instability. More specifically, we found 
that the average luminosity scales as cx B 3 with the magnetic- 
field strength, in good agreement with the expectation that the 
radiated energy should scale as cx B 2 , while the duration of 
the emission should scale as oc B^ 1 . In this way we were 
able to perform a direct comparison with the observations of 
SGR 1806-20, finding a very good agreement with the dura- 
tion of the burst and an emitted luminosity which is about an 
order of magnitude larger than the measured one. As a result, 
although our modelling is oversimplified, the overall agree- 
ment with the observations from SGR 1806-20, lends support 
to the suggestion that the basic phenomenology is that of an 
internal-field readjustment. 

Finally, we have discussed the gravitational-wave emission 
which should be expected as a result of the /-mode oscil- 
lations triggered by the instability. Also in this case, our 
calculations reproduce the expected scaling behaviour of the 
root-mean-square gravitational-wave amplitude in the limit of 
weak magnetic fields, i.e. h TSS ~ B 2 . This important re- 
sult allows us to extrapolate with confidence our estimates 
to smaller magnetic fields than those covered by our simu- 
lations and conclude that the gravitational-wave signal from 
/-mode oscillations will be undetectable for realistic values 
of the magnetic field, i.e. B p < 10 15 G and marginally de- 
tectable by third generation detectors for unrealistically large 
magnetic fields, i.e. B p > 2 x 10 16 G. 

As a self-consistent solution to this pr oblem in fully resis- 
tive r egime is close to be within reach dDvonisopolou et all 
120121) . we plan to extend the investigation reported above 
and thus remove many of the approximations that our present 
analysis had to sustain. It will then be possible to set even 
more precise connections between the complex phenomenol- 
ogy observed in magnetar flares and the violent dynamics of 
hydromagnetic instabilities in relativistic stars. 
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Figure 11. Left panel. Total magnetic energy (normalized to the initial value) versus time for simulations with B p = -Be. 5 and different values of resistivity: 
?7o/M Q = 0.02 (dashed red), 0.06 (long-dashed purple), 0.10 (blue), 0.14 (black). In the insert we have also rjo/M Q = 0.12 (dotted black). Right panel. 
Ratio of toroidal and poloidal magnetic energies versus time for the same collection of simulations. 
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APPENDIX 

THE ROLE OF RESISTIVITY IN THE ATMOSPHERE 

In this Appendix we discuss how the results of our simulations depend on the atmospheric value of resistivity 770 and how a 
suitable choice for r)o can considerably limit the influence of this free parameter. We recall that if the value of 770 is too small, we 
expect that the exterior field is not evolving rapidly enough relative to the interior dynamics, thus accumulating magnetic field 
distortions in the external layer of the star. As 770 is increased, the timescale for the exterior field evolution decreases until it 
becomes comparable or shorter than the timescale of the internal evolution. At that point the overall dynamics is less sensitive to 
a further increase of rj . Finally, if we further increase 770, it eventually will become too large and the effects of resistivity in the 
outer layer of the star will become dominant, significantly influencing also the interior dynamics. 

In Fig. QT|we show the evolution of the total magnetic energy and of the ratio of toroidal and poloidal energies, obtained 
with different values of 770 ranging from 0.02 to 0.14 M©, with an initial magnetic field strength B p = Bq.§. We first focus 
on the total magnetic energy (left panel of Fig. [TTJ and note that in the range 0.02 to 0.10 M©, the differences among different 
evolutions become smaller and smaller as the value of 770 is increased. Furthremore, while a value of 770 = 0.12 M© still gives 
comparable results, an additional increase leads again to significant differences. As a result, for these magnetic field strengths, 
770 ~ 0.06 — 0.12 M© represents a reasonable value for the resistivity. 

Let us now consider the toroidal-to-poloidal energy ratio (right panel of Fig.fTTI). which is particularly sensitive to the value of 
resistivity. When considering 770 ranging from 0.02 M© to 0.10 M©, we have a significant increase of the toroidal field production. 
Around 0.10 M©, however, the ratio becomes almost independent of 770, with no significant differences for 770 between 0.10 and 
0.12 Mq. However, if we further increase 770 above 0.12 M©, we again find a significant dependence, this time with a decrease 
in the production of toroidal-field. 

In conclusion, for 73 p = 736.5 a resistivity value of 770 = 0.10 M© yields optimal results and the scaling of this value for 
different magnetic field strengths can be done simply as to 770 /M© = 0.10 x {B p /Bq^). In the limit of low magnetic field, a 
change in B p accompanied by this rescaling of ?/o would give exactly the same dynamics, only rescaled in time. 



